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-IS,  ABSTRACT 

Using  contour  integration  techniques,  we  have  calculated  the  transient  response  of 
a  lossy  homogeneous  isotropic  plasma 'to  a  unit  step  and  a  step-carrier  sine  wave. 
We  have  found  that  in  the  presence  of  losses  the  temporal  step  function  response 
within  the  plasma  does  not  approach  zero  for  large  time  as  for  the  collisionless 
case,  but  rather  approaches  unity.  The  rate  at  whichthe  transient  approaches 
the  asymptotic  value  of  unity  is  strongly  dependent  on  the  losses,  being  quite 
rapid  for  large  losses  and  exceedingly  slow  for  very  3mall  losses.  This  result 
has  an  important  impact  on  the  propagation  of  EMP  through  the  plasma  sheath 
around  a  reentry  vehicle,  and  through  the  ionized  region  near  a  low  altitude 
nuclear  fireball. 
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Abstract 

lining  contour  mtegratiomtechniques,  we  haveealculated  the  transient  response 
of  alossy  homogeneous  Isotropic-plasma  to  a  unit  step  and  a  step-carrier  sine 
wa.e.  AVe  have  found  that  in  the  .presence  of  losses  the  temporal  step  function 
response  within  the  plasma  does  not  approach. zeros! or  large  time  as  for  the 
collisionless  case,  but  rather  approaches  unity,  'x’he  rate  at  which  the  transient 
approaches  the  asymptotic  value  of  unity  is  strongly  dependent  on  the  losses, 
being  quite  rapid  for  large  losses  and  exceedingly  slow  for  very  small  losses 
This  result  has  an  important  impact  on  the  propagation  of  EMP  through  the  plasma 
sheath  around  a  reentry  vehicle,  and  through  the  ionized  region  near  a  low  altitude 
nuclear  fireball. 
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Transient  Signal  Propagation  in  Lossy  Plasmas 


T.  INTRODUCTION 

There  lias  been  a.great  deal  of  research  donedHaskell  and  Case,  1966a,  fe, 

1967;  Knop  and  Cohn,  19G3;;Knop,  1965;  LighthiU,  1965;  Fclsen,  1969;  Antonucci, 
1972)  on  the  transient  propagation  of  electromagnetic  signals  through  lossless 
plasmas,  and.  the  results  are  quite  appropriate  for  the  calculation  of  EMP  through 
the  undisturbed  ionosphere.  For  transmissiomoFEMP  Uirough  the  piomhct  sur¬ 
rounding a  reentry  vehicle  or  through  the  ionized  region  hear  a  low  altutid-  nuclear 
fireball,  however,  the  theory  is  quite  inadequate  sinci-h  re  the  effect  of  coUisional 
absorption-on' the  transient  propagation 4s  quite  significant.  The  effect  »  collisional 
losses  lias  been  considered  .previously  by  Field  (197i?  using  the  method  of  char¬ 
acteristics,  but  unfortunately  he  does  not  present  results  for  the  transmitted  field 
iiuthe  cases  of  interesFto  us,  In  this  report,  wenviU-  therefore  extend  the  previous 
theories  to  include  the  effect  of  coUisional  absorption  on  the  propagation  of-  tran- 
sients  in  a  cold*,  homogeneous  isotropic  plasma,  V/e  shaU  consider  only- two  types 
oftransient  signals:  the  step-carrier  sine  wave  and  the  unit  step,  The  transient 


_ _ *  The  oeld  plasma. model  i3:apprpximately  valid'i3.  those  regions-near  a  nuclear 

burst  where  the  electron  temperature  is  less  than  30,  UUUCK.  Also,  in  applying 
these  results  to  EMP  propagation  v/e  must  be  sufficiently  far  from  the  firebaU  to 
neglect  nonlinear  effects, 

(Received  "or  publication  s^  April  1073} 
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response  within  the  plasma  .o  the  unit  step  is  quite  important  since  the  response  to  an 
arbitrary  EMP  can  always  be  synthesized- as  a  superposition  of  the  unit  step  responses. 
As  a  starting-point  for  cur  investigation  we  consider  the  Maxwell  equations 
for  the  electric  field  strength  E,  and  the  Langevin  equation  for  the  electron 
velocity  v.  Ifrwe  write  E  in  terms  of  the  Laplace  transform 


E(x,  t) 


E(x,p)ept  dp  , 


it  is  readily  shown  that  for  x  >  0 

x 

E(x,p)  -  &(x=0,j>)e  'c  , 


where 


’(*%) 


1/2  ,  ,  1/2 
(p  +vcP+(jp)  , 


Up  =  electromplasma  frequency, 

vc  -  electron  -  neutral  collision  frequency. 

We  have  also-assumed  that  E(x,  t=0)-0  for  all  x>0.  In  order  to  evaluate  E'x,;t), 
it  Is  necessary  to  perform  the  integration  along  the  Bromwich  contour  E  in 
E<j,  (I),  This  may  be  accomplished  either  by  integrating  dr-ectly  along  the 
Bromwich-contour,  or  by  deforming  that  contour  Into  the  left  -half  p-plane.  This 
latter  choice  leads  to  integrals  which  rapidly  converge,  and  will  therefore  be 
purs  v  d  In-this  report. 

2.  HR* 1  eiUROlNTS  ANDRRANCil  CUTS' 

The  action  y  inEq.  (3)  can  be  written  in  product  form  as 


where  s  f  ‘  '2,  '.'2  -  (p^^c)'5/2,  73  -  (p-p.>1^2,  74  s  <p-P2>1^  For  *,<-2^ 
we  have 


*  4-»2cm  1/2 , 


-i 


l 


p 

I 

%-  '  T 

I 

and  p- co,  and  we  choose  the  branch  cut  shown  inFigure  lb  for  thisquantity.  ^ 

For  yn,  whichihas  branch  points  at  p-=  and  p -=»  ,  we  choose  trie  branch  cut  1  -.f: 

illustrated  in  Figure  lc,  while  for  y^,  which  has 'branch  points  at  p  -  pg  an^  P  "  *  s-  ' 

we  choose  the -branch  cut  shown  in  Figure  la.  If  we  write  p  =  exp(iOj),  > 

(p+vc)  =  r-g  expiiO^),  (p-pj-)  -  exptiGg),  (P-P2)  *  r4  exptifl^},  we  can  write  -4 

Yl  »  ±r1l/2  exp(i01/2)J  y2  =  ±r'|/'2  exp(-i02/2),  73  *  i  r^2  exP(i03/2)  and  J 

7.  =  ±r^2exp(i04/2).  For  each-of  these  functions  we  then  choose  the  itiemann  |*- 

sheet  corresponding  to  the  upper  sign,  so  that  v 


exp 


-■(vvwy 


(5) 


and- the  branch  cuts  for  the  total,  function  7  -  are  those  shown-in 

Figure  2  for  v.Q<.  2t>0>  and  in  Figure  3  for  ^  >  2Wpi  The  integral  in-Eq,  (1)  can 
now  be  performed  by  deformingsthe  Bromwich  Contour  B  in  Figure  2  into  the 


Im(p) 


Figure  2.  Branch  Cut®  r  the  Function  7  when  v,  < 


2=  r  : 


f  ^ 
:  -s 
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Figure  3.  Branch  Cuts  for  the  Function 7  When  -v,  >  2up 

le  f t  -  haifp  -pi  an  e .  If  we  assume  that  Etx-O,  p)  has*:pole  at  so  and  Sp  thcnEq,  (1) 
can  be  rewritten- 


E(x,  t)  --  }  *  J  *  I  *  / 


E(x-0,p)=e 


where  the  contours  Cp  C2,  Cp  and  C4  are  shown  in  Figure  4. 

Toi'Justratehow  the  contour  integrals  ave  performed,  let  us  consider  the 
integral, along  Cj  in  Figure  4,  At  a  point  A  on  Cj  v<e  have  =0,  r2  *  uc~crr 
r,  *  -  (o2  -  i/co+<r2)i/2,  and  01  =  *,  02  -0,  63  +  0A_  -  0.  At  the  point  B  on  Cp 

r  r»  r*  and  r„  remain  the  same  but  is  now  -<r.  Therefore  at  point  A  we  find 
upon  using  these  results  iniEq.  (5)  that  7  *  i{o/(s'c"®)]  (up+o  '  ,  while 

7=  -i[o/(rc-(j)]l/2  (i^+e2-^}1'2.  As  a  result  the  integral  along 
contour  C^  canbe  written  as  [using  the  faet  thaf  p  -  oexpUs)  and  dp  1  -do] 

_  !p_  *  wzs“  r  d<r^(°' p  r- -a)  c_0t  expf'if(^') 1  /2  j 

-  doE{o,p=-o)e'ot  expji|,|^/  > 

.£-/*  d«£(o.  (»S2-V>>]  1,2  ■  ™ 


Figure  Deformed  Contour  of  Integration  for  the  Case  When  v,  <  2u^ 

The  integral- along  the  contour  C2  is -evaluated  in  the  same  fashion^  while  the  in- 
tegration:aIong  the  contours  C3  and  C^  is  trivial.  For  the  interested  reader,  the 
details  of  the  integration  along  C0  are  contained  in  the  Appendix  BHo  this  report. 

,3.  NPMERIC-U.  RESULTS 

3?L  Mcp-fiiimfjr  Siac-ftsiir 

*We  first  consider  the  ease  when  E'fei),  t)  is  given  by 


ftj  ■  y/2,  02  ’  tan  lujvj,  9%  -  tan  [  2(uQ-Z)/i'c)  and  <?4  =  tan  :[«2(u 0+Z)/vcJ( 
we  get 


so  that  the  integral1  over  Gg  and  C^  yields 


A  A  A  . 

°2  +  03  4  °4> 


I-  +  Ir.  =  expi(-r“F  cos  <J>  )  sin  (u i  t  -  F^-sin  <t>  )  . 
c3  t4  o  c 


Using  Eqs.  (6),  (7),  (10)  and  (A5),  we- have  calculated- the-time  response  of -the 

electric  ficld-in  the  plasma  for  the  step-carrier  sine  wave-excitation  of  Eq.  (8), 

Figure  5  shows  the^response  at -  1  for  an. overdense  plasma  (Wp/«0  "  2)  for 

different  values  of  We  comment  that  the  resists  for  vc/u  -  0, 01  are  nearly 

identical  with  those  obtained  by  HaskelFand  Case  (1966a)  for  the  case  J^/w  -  0. 

For  v  jtj  =1.  however,  the  response  is-quite  different.  This«difference  is  even 
C  P  (JaX  , 

more  apparent  for  the  cases  when  — -  -  5a»d-^—  -  10,  whichiareshown  in 

Figures  6  and  7.  For  large  values  of  t^x/.c,  even  when  rc/w  -  0.1,  the  transient 


1 

(■m 

[tJEJ 

mm 

Figure  5.  Transient  Response  to  Step-Carrier  Sine  Wave  for  «  x/e  -  1  and 


response  is  alteredsignificarllv,  and  for  vc/u^  -  1  the  character  of  the  response 
is  completely  different  from  the  collisionless  limit.  Figures  8  and‘9  illustrate 
howthenature  of  the  transmitted  signal  changes  as  the  distance  into  the  plasma  is 
varied;  We  note  from  Figures  8  and  9  that  for  large  iQx/c ,  there  is  an  oscillation 
followed  by  a  long  wake  in  which  the  amplitude  of  the  signal  isnearly  constant. 

It  is  readily  shown  that  in  this  wake  region,  E(x,  t)  can  be  approximated*  by  (see 
Appendix  B) 

E(x,  t)  exp  (-F^-cos  <h  )  sin  (uQt  -r-^-sin$)  + 


In  deriving  Eq.  (11)  we  have  assumed;  that  t  »  x/c  and  that  v  /u  is  of  order 

c  p 

unity.  From  Eq.  (11)  we  see  that  Tor  largel  the  transient  signal  consists  of  a 
sinusoidal  wave  at  frequency  . u  plus  a  component  whic-ft'eventually  decays  tor/ards 
zero;aS‘t 

Figures  lO  through  12  indicate  the-natiire  of  the  response  for  up/uo  s  4‘, 

Figure  10  shows  the  transient  response  for  uQx/c  -  0.  5,  where  the  effect  of  large 
values  of  t>c  is  »o  radically  alter  the  signal  amplitude.  From  Figure  11  we  see  that 
when  uDx/c  =  2.5,  both  the  amplitude  and  the  character  of  the  transient  response 
are  radically  changed  from  the  case  of  zero  collision  frequency.  In  fact,  the 
frequency-of  the  oscillation  for  vc/u  -  1  is-quite  drastically  reduced  from  the 
case  when  i>c/up  ~  0.01.  Figure  12  illustrates  the  transient  response  at  u0x/c  -  5, 
Atethis  distance,  for  -  1,  there  is  hardly  any  oscillation  of  the  transient 

response,  but  rather  a  slowbuildup,  and  then  an  eventual, decay  (not  shown  in  the 
figure)- in  accordance  with  the  second  term  inJ2q,  (Id),  since  the  first  term  in 
Eq,  (11)  is  quite  small  for  this  case. 

In  Figure  13,  we  have  shown  the  effect  of  collisions  on  propagation  in  an 
underdensa  plasma.  Here  we  note  that  no  interesting  or  unexpected  behavior 
occurs,  and  we  wilt  not  study  the  underdense  case  further?  In  addition,  this  case 
has  becn  previously  considered!)}'  Case  and"  Haskell  (1967).  We  only  comment 
that  our  results  do  agroe  with  the  previous  ones  obtained'by  Haskell  and  Case 
CiB66SMn;the  ilmtt  t'c/up  -s  o. 

';To  check  the  accuracy  of  Eq.  (11),  let  us  consider  the  case  when  i'  /w  =  l, 

c  p 

w  [u^  -  2  and L/c  -  10.  If  we  take  a  point  in  the  signal  wake  at  u^t  -  15,  we 

hive  from  Eq.  (U)  that  E  -  2.  54  X  10*3,  while  the  exact  calculation  (see  Figure  9) 
gives  E=  £.51  X  10‘3. 


3;2  Step  Function 


We  next  consider  the  situation  when  K(x -0,  t)  is  given  by 


E(x=0,t)  --  u(  ti 


(12) 


For  this  case  the  pole  sQ  coincide**  with  the  branch  point  at  p  -  0,  while  there  is 
no  pole  Sj.  Here,  special  care  must  be  takeit  in  performing  the  integral  Iq,,  and 
the  integral  in  Eq.  (7)  becomes 


(13) 


where  I.  is  the  integral  along  a  circle  of  radius  6  surrounding  the  origin  in  the 
p -plane.  That  is 


f  d0e6t  exP(i0)  e 


7(P -6  exp(ifl)) 


(14) 


Using  the  result  of  Eq.  tl3),  plus  the  appropriate  expression  for  Ic  from 
Appendix  A,  we  have  calculated 

f  f  a  (pt  “7) 

E(x,  t)  *  J  +  J  E-(x=0,.p)  e  c  dp  .  (15) 

Ci  C2 

Figures  14  and  15  show  the  transient  responses  to  the  unit  step  input  at  distances 
UpX/c  “  l  and  a  within  the  plasma  for  rc/up  -  0.  01,  0,1,  and  1,0,  It  is  interesting 
andimportant  'o  note,  from  Figure  15,  that when  v/u  =  0. 01  the  transient 
response  for  qpt  <  40: is  nearly  identical  to  the  results  obtained  by  Haskell  and 
Case  dOSGaWor  vQ/Up  -  0,  However,  we  now  note  that  when  t>c/wp  /  0- the 
transient  response  doesmct  approach  zero  as  t  —  co  as  it  does  for  rc/u i  -  0,  but 
rather  approaches  unity.  It  is  clear  from  Figures  14  and  15  that  the  rate  at 
which. the  transient  grows  lowards  unity  js  highly  dependent  on  the  collision  fre- 
quencyi  For  example,  if  v  fu  -  1  we  see  frorn  Figure  15  that  the  Iranniont  has 

*  up 

reached  a  value  of-e  by  upt— 16;  for  vc/up  ■  0.1  the  e  value  occurs  at 

60;  while  for  v  /o  *  0,  01  the  signal  has  reached  of  its  final  value  of 
E*  -  v  P 

unity  at  wpt—  1 ' ?05>.  Figure  16  illustrates  the  transient  *  esponse  at  the  Larger 
distance  tppx/c  *  10,  The  character  oithe  response  is  essentially  the  same  as 
for  u  x/c  «  5,  exccpt  that  the  rate  of  growth  is  toward  the  asymptotic  value  and 
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Is  even  slower  than  for  u^x/c  =  5,  T  s  point  is  furth'r  illustrated  in  Figure  17, 
where  we  compare  the  responses  for  t,  x/c  =  5,  u  x/c  =  10  and  u  x/c  *  20,  for 

F  sf  r 

the  case  when  v_/w  =  JL 
c'  p 

The  behavior  exhibitedinFigures  14  through  17  is  j,  st  what  is  expected 
physically  since  a  wave  of  frequency  w  propagates  as 


T  ,  2 

,  i  fl  -  -pj 
c  l  u(u+ 


The  behavior  near  the  leading  edge  of  the  responses  ’s  determined  by  the  high 
frequency  components  whfch  propagate  as  exp  (i-~~),  while  thelime  asymptotic 
behavior  is  determihed  by  the  low  frequency  components,  which. propagate  as 

f  /  \  /,  ? , .  v 


which  approaches  unity  as*w  o  . 

It  is  possible  to  develoi  ail  accurate  analytical  expression  for  the  field  strength, 

which  predicts  the  time  asymptolic  behavior  we  have  observed  in  Figures 

14-through  17,  Lot  us  assumethaf  t  »  x/.c  b  t »  1,  and  r  /w  is  of  order  unity,. 

P  C  P 

Then  it  is  shown  in  Appendix  C  that 


E(x,  t)  -  1  -  Erf 


2c(vct)J 


This  expression  is  quite  accurate-for  t  »  x/c  and  vc/Vp  ~  0(1).  To  illustrate  this 

we  have  compared  the  results-of  Eq,  (16)  withlhe  exact  result  in  Table  1,  for 

v  jw p  =  1  and  UpX/c  =  1  and  ID,  Wenote-tliat  the  agreement  is  excellent  once  we 

are  not  too  close  to  t  =  x  (e.  EquationftGVis  not  nearly  so  accurate  for  smaller 

values  of  vJu  since  then  it  lsmolonger  legal  to  neglect  a  in  comparison  with 

v  in  the  argument  of  the  sinefunctionlq  Eq,.  (13),  A  comparison  of  the  exact 

results  with  the  approxin-ation  of  Eq,  (16)  for  v  jo  -  0.1  is  given  in  T  .oie  2, 

c  p 

As  expected,  the  agreement  ismot. generally  good. 


i.  DISCUSSION 

The  most  interesting  point  we  have  noted  in  tbls  report  Is  the  crfect  of  plasma 
losses  on  the  transient  response  of  a  unit  step.  As  pointed  out  in  Section  1,  this 
result  has  Important  consequences  for  the  interaction  of  EMP  with  the  plaama  . 
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surrounding  a  reentry  vehicle*.  For  example,  suppose  the  plasma  sheath  thick¬ 
ness  x  =  10  cm,  u  =  1.  5  X  to10  sec"1  and  r  =  1. 5  X  1010  sec““.  Then  v  fu  -  1 
p  c  c'  p 

and  WpX/c  =  5.  Now  suppose  we  had  a  square  electromagnetic  pulse  or  duration 

TQ  =  10~®  sec  incident  upon  the  plasma  sheath  surrounding  s  reentry  vehicle. 

Then  from  Figure  16  it  is  clear  that  the  nature  of  the  signah  present  at  the  skin  of 

the  reentry  vehicle  when  is  significantly  different  than  for  r  =0.  This 

difference  is  illuscrated  qualitatively  in  Figure  18.  It  is  possible  using  the  results 

we  have  obtaihedtfor  the  unit  s-tp  exciiatioh-to  calculate  the  transient  response  for 

an  actual  EMP  due  to  a  nuclear  blast,  but  in.  order  to  keep  this  report  unclassified 

we  have  not  presented  this  result  here, 

The-calculatiom  of  the  transient  response  presertced  here  is  i  eadily  extended 

to  spatially  inhomogeneous  plasmas,  by  approximating'the  plasma  bv  a  series  of 

homogeneous  layers,  each  with  different  value s  of  u  and  v  .  Then  in  the  nth 
A  P  c 

layer,  E(x,  p)  has  the  form 

X 

En  =  Ah(p)>eTnc  +  Bn(pleYnc  (17) 

where  An  and  B(1  are-obtained  by  requiring  that  E  and  3E/Sx  be  continuous  across 

the  interfaceibetween  each  layer,  and  y  is  given  by  Eq,  (3)  with  the  values  of  a 

n  p 

and  i>c  in  the  nth  layer  used. 

The  methods  used  here  can  also  be  readLly-extendedcto  calculate  transient 
behavior  in  a  lossy  dielectric  or  in  a  conductor.  For  this  case,  it  is  readily 
shown  thatnhe-stepTunction  response  at  a  point  x  within  the  material  is  given  by 


E(x, t)  -  l  -  —  f  e"^rsin  iy^U-y)1/2  , 

'  Jo  y  L 


where  r  =  (o0/c)t,  i  =a0x/(ev),  v  =  (pQc)-1^2,  e  is=the  dielectric  permittivity 
of  the  material  and  aQ  it's  conductivity.  For  t  »  1  we  obtain  as  an  asymptotic 
representatipsiifor  the  step  response 


E(x,  t)  -  I  -  Erf 


so  that  ir.  a  cond  Xvor  or  lossy  dielectric,  the  unit  stop  function  responsc  at  a 
given  point  in- the  material  eventually  approaches  unity  for  large  values  of  t.  A 
plot  of  tLe  transient  respcn=;  to  a  unit  step  in-a  lossy  dielectric  or  conductor,  as 
cornputsd  from  Eq.  18,  is  shown  in  Figure  19, 

*It  is  of  little-consequence  for  the  EMP  propagation^  the  ionosphere  since 
there  is  so  small  that  v,/u i  «  1, 


Figure  19.  Transient  Response  of  a  Conductor  of  Lossy  Dielectric  to 
a  Unit  Step 
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Appendix  A 


Here  we  discuss  the  integration  along-the  contour  C9.  Consider  the  point  C 
shown  in  Figure  3.  Here  tj  5  Pl  -  ii^+v^+Jh1^,  rg  S  p2  =  (n2-vQrj+u2H2  , 

r3  E  -63  =  v'  r4s  ^4  1  (4up-1/c+r,  )  and  °l  s  0l  - /2  *  ‘an-1  [U'c/2+rj)/ZJ» 

°2  s  ®2  ^  tan_1 1  ZHVJ2  'n)J  «  ®3  Si%  =  04=  T4.  s/2  +  tan-1  (n/2Z).  At  the 

point  D,  all  the  above  quantities  are  the  same  except  03  which  now  equals  -s. 
Defining  6>c  *  1/2  and  E  -  <p1p3p4/p2)^“,  we  can  then  write  the 

contour  integral  on  the  upper  horizontal  portion  of  C2  as 


1  j  -(n+"^-iZ)t  'i'tsinOc)  A 


E(p=-w=-n+iZ) 


y 


,  ~R{cosO  4.£sin0  ) 

+  1“  dne  e  c  c 


Ete^-^-jrHZ). 
*  «: 


The  integrals  over  the  lower  horizontal  portion  of  C2  ca.rbe  performed  in  a 
similar  fashion,  and  yield 


<>  i  ~<n+^+iZ)t  ^R(ccsfl  -isinfl) 

h‘-m  yo  'd^e  [p 


-~R(cos  6  ■ 
c  c 


isine)1  A  i' 

c  ;  E(p  -  -~f  ~  n  -  iz)  . 


Combining  Eqs.  (Al)  and  (A2)  we  may  write 


■  «  C-—  f2  +!  -nt  XRcos  °c 

h  +  h--  *  ;/0  »  °Im 


(a  V 

<  E(p=-|-i?+iZ)e 


i(Zt-|RsinGc) 


V 

—£4.  F 


2 4,5  -nt1rRcos*c 

^  dn  e  ”  e  CX 


J2  d* 
Jo 


(  a  J7..  i(Zt+£  R  sin  0  )) 

j  Elp=~^-n-iZ)e  c  CS.  (A3) 


Finally  we  must  perform  the  integrals. over  the  vertical  portions  of  C2.  Consider 
the  point  G.  At  this  point  a  3^  =  [(^+5  )2  +-p2]1^2,  r2  =  r2  =  (f^p2)1^2. 


r3  s  V  [(‘'Jt+V2  +  (Z-ft)2  11/2.  r4  «  5t4  =[{vj2  +  tf  +  (Z+p)"j1^  , 

%  "  %s  r  -tan'1(.p/(vc+5)),  p2  -  ©2  s  ? -tan‘rl{p/5 ),  03  -0^=  -s+tan-1 
{ (Z-p)/(  -f-+5  )].  04  =  04  '  *  -  tanr-[;(Z^)/(^4-5)j. 

At  the  point  II,  all  the  above  quantities  remain:  the  same  except 
03  *  v  +  lan‘1[<Z-p)/(*jr+  5)]  •  Defining  *-  +04)  and  R0  = 

(r,r3r4/r2)^2,  the  integrals  on  the  vertical-portion  of  C2  can  be  written  as 


2,1/2 


pn  _  e 
2  ‘  “ 


■<  c+?)t  F  cos*  I A 

_ _ -  f  fine  c  0  0  R 


]  dpe 

•0 


Rep(ps“V5+ip)e 


i(pt-|RoSin^)J 


Z  cos  0 

f  dp  ec  0  0  Re 

J0 


(a  Hpt  + 

7  E(p--rc-5+jp)  e 


rRosinV 


The  total  integrat  on  the  contour  Cg  is-then  given  by 


f  ft  It* 

L,  *  I*  *  I  + 


Appendix  B 


Here  we  outline  the  derivation  of  Eq.  (11).  We  consider  the  limit  when 

t  »— .  u  t»  1  fid  f  !u  is  of  order  unity.-  In  this  case  It  can  be  shown  that 
c  p  c  p 

the  integral  alonj;  .  he  contour  Cg  is  negligible  so  that 


E(x,t)  I„  +  lr  Up  , 
n  *"3  S 


wherein  isgiven  in  Eq.  (7)  and  I,.  +1,.  is  given  in.Eq,  (10),  Now  in  Eq,  (7) 

C1  ^3  '“4 

it  is  clear  that  for  t  large,  the  principal  contribution  to  the  integral  must  come 

from  a  near  zero.  Therefore  I,-,  may  be  approximated  by  expanding  the  integrand 

a 

in  Taylor  series  aboute  -  0.  We  get 


1  /•* 

U  ct  -i-  f  da 
e  i  ■"'“o  v 


ix) 


1/2  upx 


w0  •'q 


f  doe" 


iin  J2_, JJUi 
V"c  / 


The  integral  inJEq.  (B2)  Is  a  standard  Laplace  Transform  and  yields 


1  Zuc(Vr) 


A  e  v  t 


Appendix  C 


For  a  unit  stop  when  "c/wp  is  of  order  unity  and  t»  x/c  we  may  neglect  Ic2 
so  that  the  transient  response  is  given  by  Icp  which  from  Eq.  (13)  with  5  0 


E(x,  t)  *  l  ~  -^e'uT  sin  L  ^H-y2-Xy^  ^  ,  (Cl) 

where  T  -  u  t,  X  -  rc/up,  L  =  upx/c.  When  T»  1,  it  is  clear  Shat  the  principal 
contribution  to  the  integral  comes  from  y  near  zero.  Expanding  the  integrand 
about  y  =  0  ana  extending  the  range  of  integration  to  »  tlien  yields 

“fay'*  ■  ,c2> 

The  Integral  in  Eqt  (B2)  is  a  standard  Laplace  transform  (Erdelyi  et  a’,  1954) 
and' upon  evaluating  it,  we  gel 


E(x,  t)  -  1  -  Erf 


„2{XT)1/2  ] 


Preceding  page  blank 


''  '  '  .  -  Appendix  D 

Here  we  include  a  Fortran  listing  of  the  computer  program-used  to  calculate 

thetransient  response  of  a  lossy  homogeneous  plasma.  The  program  inputs  are 

XL  =  wpx/c,  AL  -  vju p,  OM  *  w0/up.  D  =  the  radius  5  of  the  contour  in 

Eqs>  (13)  and  (14)  (D  is  usually  set  to  0.  01),  M  -  number  of  different  values  of 

T  =  u  t  at  which  calculation  is  to  be  carried  out,  and  L  =  1  for  a  unit  step  input, 

P 

while  L  =  2  for  a  step-carrier  sine  wave  input.  The  program  outputs  are 
E  E  =  E(x,  t)  =  electric  field  strength,  F  F  -  contribution  to  E  from  poles, 

FU  -  contribution  to  E  from  integral;  along  C^,  and  FUNC  -  FUN  =  contribution 
to  :E  from  integral  along  Cg. 
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PROGRAM  RON  ( IN°UT ,  OUTPUT;) 

REAL  M2-,  42  ,  H4,  Nl N£,  N3,  N4 
COMMON  AL,  OH,  XL,  l,  PI,,  ZO 

COHMON/TERMS/TERMl,  TERMiSI,  TERMi ,  TERM2SQ,  TERH3 ,  TERM3.3Q 
COMMON/3/9 

5  READ  13,  0,  AL,  OH,  XL,  L,  M 

IF(AL  .EQ.  U . 0 j  STOP 

lu  FORMAT f  4 Flu.  4,2:14) 

PRINIT  2u,  AL,i  OM,  XL,  L 

2u  FORMAT  (1H1,49X‘»  *AL  =  *Fi0.4/50X, *0M  =  *lrlU.4/50X,  *XL  =  *c10. 
1,*  v  =  »IS//-7SX,*f*,l3X>=*cE*,ioX,*FF*,16X,»-0*,l5X,»FUN»,15> 
2*FUNCS  7X,*CP  3PCCNQS*/) 

PI  =  3. 14159265353079 

TERHi  =  9.u  *  AL 

TERH1S3  =  TERHI  *-*2 

TERM t  s  10.0  ’•  AL 

TERN2SX  =  TERH2**2 

TERM 3  =  9.5  *  AL 

TERM3S1  =  TERM3**2 

ZO  =  SRRTtl.O  -  j;25*AL-**2> 

Pd  =  3HT(AL**2  +■  OH* *2) 

M2  =  AL  /  F2 

F3  S  SaRn  (OH  -  ZO)"**^  4  (AL  /  2. 03  **2) 

M3  =  j.5*AL  /  F3 

r 4  =  SRRT (  (OH  t  Z0)**2  1  (AL  /  2.0)**2) 

M4  =■  d.5*AL  /  F4 

G  =  S^RT  (OH)  *  SIRT (FJ)  *-  33RT(c4)  f  S0RT(F?) 

PI  =  l.u  /  SORT (2.0) 

HI  =  l.u  /  SQRT(e.O) 

P2  =  3RRT ( U . 5*( 1  -  M2)  ) 

M2  =  5IRT  (  0.5*(i  +  H 2)  ) 

M3  =  5IRT(J.5*(1  f  M3)  ) 

P3  =  STRT (0.5*( 1  -  83)  ) 

IFCOM  . LT.  ZO)  P3  =  5°3 
M4  =  SRRT ( 0.5*( 1  f  rf4)  ) 

P4  =  3XRT(0.5*(1  -  H4)  > 

CP  =  N1*N2*N3*N4  *  Ni*.N4*P2*P3  -  N2*M4*P1*P3  *  N3*N4*P1*:>2 

1  -  N2*N3*P1*°4  *■  Ni*M3*P2*=4  -  N1*N2*P3*P4  -  P1*P2*P3*®4 

S P  =  Pl*N2*N3*N4  -  P2*Nl*N3‘M4  P3*M1*N2*M4  f  P1*P2*P3*N4 

1  f  ®4*Nt*M2*H3  >  N1*P2*P3*?4  -  N2*P1*P3*P4  *  M3*®1*P2*®4 

DO  30  J=i,M 


Quoigg 

oobiib 

000120 

000130 


000150 


XX  s  FLOAT (J> 

IF (XX  ,lE.  6.0)  T  =  XL  +  J.02*(XX  -  1.0) 
rc(XX  ,:3E.  7,0  .  AMO.  XX  .LE.  21.  J)  T  =  XL  *  0. 
IF(XX  ,5E.  22.u)  T  =  XL  +-3.1  *  (XX  -  21.0) 
IFfL  .EO.  1)  FF  =  FJMR(T) 

IF (L  .E3.  2)  FF  =  BXP(=-G*X.*CP)  *  SIN(QM*T  -  3 
FUNFU  =  FU(T) 

FUNFUN  -  FON  (T) 

FUNFUN3  =  FUNC(T) 

EE  f  F*  f  FUKFSI  --  F4JNFUH  t  “  JMEUNC 
3PT-TMS  =  SECOND (A) 

30  PRINT  40,  T,  EE,  FF,  FUMF-J,  FUNFUN,  FJN=DNG,  3 
40  FORMAT ( 5X, r 8 . 3, IP, 5(3X,E15.3) , OP ,3X ,F9. 3) 

GO  TO  5 
EMO 

*¥*,.¥**,.**¥**¥*»**¥*-*^*  ********  ********  **,************ 


000170 

X ,  *XL  =  *c10  •  4/5uX  000=190 
15X,*FUN*,15X,  00019C 

Q0Q2uu 
000210 
000220 
000230 
U00240 
00u25C 
000250 
000  270 
0QU280 
000290 
GOOlUO 
0.,0210 
000320 
000330 
000340 
000350 

ooMso 

Dud  37  0 
000380 
000 39G 
000400 
000410 
0  0.0  420 
000,430 

M3*N4*P1*32  OuO 440 

P1*P2*P3*°4  000450 

P1*P2*P3*N4  00u460 

N3*°l*°2‘®4  bu  0-470 

000480 
000490 
000.su  6 

1  +  0.2* (XX  -  6.0)  000510 


*XL*SP) 


°TIH£ 


000520 
0.00  530 
buU-540 

000570 

0siu580 

000590 


00OF20 

000630 

M  *  V**  *************05jUg4{) 
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****** ****^.***************** ************ *****************************'MM‘00d650 


FUN3TX0N  FU<T) 

COMMON  AL,  OM,  XL,  L,  PIj  10 

COMMON/ TTiVTTi 

COMMON/O/D 

TT1  =  r 

N  =  1 

TWO  =  J.O 

GO  TO  (5  s  6)  ,  L 

5  M  -  (U.99J-  0)*AL  /  2.0 
FOUR  =  FUNKKO'AL  +•  X) 

ENOS  =  FUNK1(0*AL)  4  FUNKKO  499*  AL-) 

GO  TO  9 

6  H  =  0,W9*AL  /  2.0 
FOUR  =  FUNKl 00 

ENOS  =  F U NK1 ( 0 . u 0 U « t >  4  FUM<1  (0..199**v> 

9  SUMO  =  (ENOS  4  4.0*F0dR)  *4/3.6 
10  H  =  H  /  2.6 
N  =  2  *  N 
TWO  =  THO  4  FOUR 
FOUR  =0.0 

Y  =  H 

IF (L  .EO.  1)  V  =  0*AL  +  H 

I  =  0 

20  I  =  I  4  1 

FOUR  =  FOUR  4  F-UNK1K-Y > 

Y  =  Y  4  H  4  K 

IF (I  iLr,  N)  30  TO  2u 

(>U  =  (ENOS  t  2.0*THO  4  4.0*rOUR>  *  H  /  3.0 
IF(ABS(SUMO  -  FU)  , L‘T .  1.0E-6)  RETURN 
SUMO  *  FU 
GO  TO  10 
END 


Y  =  0*1L  4  H 


000660 
000670 
0 00 690 

80U690 

0007u0. 

000710 

000720 


000760 

000770 

000780 

000790 

000800 

0j0810 

000820 

000830 

000840 

0u0850 

000870 
000880 
000890 
000900 
000910 
000920 
000930 
000940 
Di.0950 
Go 0960 


************************************************************************ 00u970 

********************************************************-»**  ************* ggogeg 


FUNCTION  FUN (T) 

COMMON  AL ,  OM,  XL,  L,  PI,  ZO 

C0MH0M/TT2/TT2 

TT2  =  T 

H  =  9.5*AL  /  2.0 
N  ■  =  L 
THO  =  0.0 
FOUR  =  rUNKR(H) 

ENDS  =  FUNK2  (0. U J  4  FUNK2(9.5*AL) 

SUHO  =  (ENOS  4  4.d*F0UR)  *  H  /  3.0 
10  H  =  H  /  X.O 

N  =  2  *  W 

THO  =  THO  4  FOUR 
FOUR  =0.0 

Y  =  H 

X  =  0 

20  I  =  1  4  1 

FOUR  =  FOUR  1  FUNK2(Y) 

Y  =  Y  4  H  4  H 

IF  (t  .LT.  NJ  30  TO  20 

FUN  =■  (ENOS  4  2 . 0*THO  4  4,a*F0UR>  *  H  /  3 
IF( ABS (SUMO  -  FUN)  .LT.  l.QE-6)  RETURN 
SUMO  =  PUN 
GO  TO  10 
END 

******************************* ***  ************* 


000990 

ooiooo 

Oul 010 
001020 
Gulo30 
001040 
UulO50 
Owlo&O 
001070 

aoiuap 

001090 
Oul TOO 
001110 
001120 
001130 
Oul 140 
001150 
001160 
Oul 170 
001180 

.0  U 01 190 

001200 

001210 

6*1220 

601230 

***********•************:*  240 
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**»*.*******♦******  ***>*****.******•*  *** <,***-*  **+*♦»***■*  *******  ******»**.****  *Uirl-250 


FUNCTION  FUNC(T) 

COUPON  AL,  OM>  XI,  L»  PI,  20 

COMMON/ TT3/TT3 

TT3  =  T 

H  =  ZO  /  2.u 

N  =  1 

TWO  =  0.0 

FOUR  =  FUNK3 (H) 

ENDS  =  FUNK3 (0*  0)  +  FUNK3(ZD) 

SUHO  =  (ENOS  +  4.0*FOUR)  M  H  /  3.Q 
10  H  =  H  A  2.0 

N  =  2  *  N 

TWO  =  THO  F  FOUR 
FOUR  =0.0 

Y  =  H 

1  =  0 

20  I  =  I  *  1 

FOUR  =  FOUR  +  FUMK3( Y) 

Y  =  Y  F  H  F  H 

-I PM  .LT.  N)  GO  TO  20 

FUNC  =  (ENOS  4-  2. 0*TM0  *■  A. 3* FOUR)  *  H  /  3.0 
IF(AB3(SUMO  -  FUNC)  .LT.  1.0E-6)  RETURN 
SUMO  =  CUNC 
GO  TO  10 


00126b 

001270 
0012/30 
0ul?9Q 
0U1-30G 
0-1310 
00 13 20 
001330 
001 340 
uo 13 50 
001360 
C  .1370 
001380 
001390 
0ul4u  0 
0-1410 
001420 
U01430 
0O1440 
0  u 1450 
UL1460 
001470 
001480 
0ul49P 


ENO  ’  001500 

*«***>  **»***.**************.******¥********»***%*:***  ************  •“>**f*****0G1510 

***»*¥  ****#*♦*.♦♦***»*  ¥.*¥*»'*♦***♦*  ***»*********.¥-*****  *»*»***<»»•►  *»*»-*->*  *»*0ul  520 

FUNCTION  FUNKi'(Y) 

COMMON  AL,  OH,  XL,  L,  PI,  ZO  001500 

COMMON/ IT  JL/-TT-1  t  001550 

tcdm  =  ?TU,YI  t  Y**2  -  AL*Y>  /  ( AL  -  Y>  )  )  Gul560 


TERM  =  5IN(XL*SQRTi<Y*(l.J  «•  Y»*2  -  AL*Y)  /  (AL  -  Y>  )  ) 
GO  TO  (10,  20),  L 

10  FUNK1  =  -EXP(-Ttl*Y)  *  TERM  /  (»I  ‘  Y) 

RETURN  „  „ 

20  FUHK1  =  OM  f  EXP(-TT1*Y>  *  TERM  /  (PI  *  (Y*»2  «■  0M**2>  ) 
RETURN 


U01570 

0olS80 

601590 

0U16U0 

001610 

U61620 


*♦♦*******.*♦*♦**¥*.¥**  *******  *****<m,  **»♦***¥***¥**¥**♦****♦********  061630 

***********************************  *********  e  ******  **¥***>****«**,**0{,},6£,g 


FUNCTION  FUNK2(Y) 

COMMON  AL,  OM,  XL,  L,  PI,  ZO 
C0H.M0N/TT2/TT2 

AA  =  SQRT( SQRT <4.0  -  Al»*2  f  Y*'2)  > 

AB  =  S1RT  (SORT (9**2  F  AL»Y  *  l.i>  ) 

AC  =  SCRT(SQRT(  Y>*2  -  AL*Y  4-  1.0)  ) 

R  -  SQRT(Y)  *  A  A  *  AB  /  AC 

AO  =  ZO  /  SORT ( Y**2  -F  l.Q  f  AL*Y) 

BO  =  ( A-_  /  2.6  r  Y)  /  SORT( Y  * *2  F  l.J  -  AL*Y) 

CO  =  2.0*20  /  SORT (Y**2  *  4.0  -  AL*“2) 

Cl  =  SRRTi  (  X. u  r  Au  ;  /  2*0) 

C2  =  SORT  (  (  l.  U  4-  BO  )  /  2.0) 

C3  =  S3RTC  (  1.0  +  CO  )  /  2. 0) 

51  =  S1RT (  (  1.0  -  AO  >  /  2.0) 

52  =  3RRT(  (  1.0  -  BO  )  /  2.0) 

53  =  SORT (  (  i. 0  -  CO  )  /  2.0) 

CT  =  -(C1FC2FC3  f  Cl*32*So  -  C2*Sl»S3  *■  C3*S1*»2) 

ST  =  -( S1*C2*C3  -  S2*Ci*Co  *  93*01*02  4-  Si*S2*93) 


001650 
001660 
0ul67  0 
001680 
0ul69C 
601700 
001710 
001720 
0  01-73'' 
001740 
031750 
001760 
001770 
001780 
0ui-79P 
001800 
001810 
001820 


„'-T  _  ^  ' 


1  tt  ^  ^tVCK  .  '{r 


S3  3  =  3XN ( ZO+TT  2 
CC3  =  33SOZO*TT2 
SS6  =  3IN(ZQ*TT2 
CC4  =  C jS(ZO*TT  2 

551  =  SiN(Z0*TT2 
CC1  =  30S(Z0*TT2 

552  =  3IN(Z0*TT2 
CC2  =  D0S'(Z0*TT2 
EXP  1  =  EXP(-XL*R*Cr) 
EXP2  =  EXP  (  XL*R*nt3 


XL*R*ST 

XL*R*3T 

XL*R*ST 

XL*R*ST 

XL*R*ST 

XL*R*ST 

XL»R*ST 

Xl*R*ST 


PI) 

PI) 

°i; 

Pi) 

PI/ 2. 0 ) 
PI/2.0) 
PI/2.0) 

pi/2-.  a ) 


Y*TT2) 


AL*Y  *  Y**2)  )  *  EX°3  * 

-  EXP2KZ0*SS2  -  TER -11*002) 


) 


EXP3  =  c  X  P  ( ■*.  AL*  TT2/&.  o 
TERH1  =  AL/2.0  ¥  Y 
GO  TO  (10,  20),  L 
10  FUNK2  =1.0/  (PI  *  M.O  t 
K£XPlKZO*SSi  -  TERHL*CC1> 

RETURN 

>V  TERM 2  =  TERM!** 2 

TERM3  =  (ZO  -  OH) *  *2 
TERM  4  =  (70  +  0Mi**2 

FUNX2  =  OM  /  (?  I*  (TE.RM2  +  TERM4)  *  { T-ERM2  +  TERM 3)  )  *EXP3*  (EXP** 
(r(Z3**2  -  0H**2  -  TERM2)  'SS3  -  2  .0*ZA*TrRMis2C3)  -  EX 2 2  * 

2*  ((ZO**2  -  OM* *2  -  tERM2)*R§4  -  2 .  u*’A*3'ERMi*0C4)  ) 

RETURN 

END 

****;*,l*«i*  **♦*¥  *¥  ♦»-♦*»***  *«**.**’****  **  **»*♦*¥  ¥*¥*+** ******  «.**♦*¥  *  *  ***** 
*  *****  ********«»)•  *  is.****  *****************  **************  ********** 

FUN3T:I3S  FUNKS (Y) 

REAL  <1,  K2*  <3,  <4 
COMMON  AL,  OM,  XL,  L,  PI,  ZO 
C0MM0N/TT3/TT? 

OOMMON/TERMS/TERMl, 

C  TER.M1  =  9.0  *  AL 

C  TERM  1-SQ  =  TERM  1**2 

C  TSRM2  =  10.0  *  AL 

0  TERM2SA  =  TERM 2**  2 

C  IERM3  ~  9.5  *•  AL 

C  TERM3SA  =  T£RHo**2  4 

YSQ  =  Y*  *2 

BA  =  3RRT(SQRTI(iERM23a  4-  Y50> 

89  =  S3RT (S3RT(TERH3S3  4-  <79  - 
0C  =  SORT (SORT (TERM3SA  +  (Z3  * 

R  =  94  *  33  *  3C 


TERMIS3,  TERMS;,  TRH2SQ,  TFRM3,  TERM3SQ 


/  (TERM1SA  * 

Y)**2)  ) 

Y) **?)  ) 


YSi)  >  ) 


A2 

32 

C2 

02 

Ki 

K2 

K3 

X4 

ft 

T2 

13 

T-4 

CS 

1 

ss 


TERH2 
TER  Ml 
TSRH3 
TERM3 
SARK 
SARK 
SARK 
SARK 
saR  T  t 
SARK 
SARK 
SARK 


SART (TERM2SQ  ♦  YSQ) 
SORT ( T£ RMiSA  ♦  YSA> 
SARKTERMiSA  *  (20 
SRPT  (TERM3SQ  4-  (20 


/ 

/ 

/ 

/ 

(1.0 

(l.o 

(i.U 

(1.0 

csi.o 

M.U 

cl, a 

( i.u 


-  Y  )*J,2A 
4-  Y ) **2) 


A2) 

32) 

02) 

02) 

mu 

32) 

C2) 

02) 


2,0) 
2.0) 
2-.0>- 

2,0) 

*  -  Mt 

/  £.UI 
/  2.0) 
r  2,3) 
/£.  0) 


PA  = 


<1*K2*K3*K4  -  <1 *K4fT 2*T3 

-  <3*<2*T1*T 4  f  <l*K3*T2*14 

-  Tl*X2*K3,t<4  ¥  T2*>X1*'<3*X4 

-  <1*K2*K3*'T4  *  <i*T£*T 3*T4 
IX p (-TERMS  *  TT3)  /PI 


K2*<»4*Tl*' 


+ 

4  K1 *<2*13*3 4 
¥  f3*«L*X2*K4 
-  K2*Tl*T3*T4 


j  +  K3*K4*Tt*T2 
+  T1*T2*T?*T4 
t  T  1>T  2*'r3*l<4 
-  K3’,T1^T2*T4 
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U01830 

Oul'OAO 

091850 

Or>1360 

001870 
ujiaeo 
6*-l  89  0 
ooisoo 
,M1°1 0 
001020 
001930 
001940 
001950 
uui-960 
0^1971 
ujlAftO 
0  0  i-99  0 
Q020uG 
0  0  2  u  1 C 
6 0202 u 
002030 
002040 
002050 
6j2v60 
*>*u02G7o 
***0j2o80 
0.42*90 

Du21(i0 
002110 
OU212Q 
Gt2i30 
0-2 140 
012150 
002160 
lo21? u 

0U2i  ef 

01-2190 
220  0 
■»?in 

41--20 
U02230 
0G224A 
302250 
0  -2260 
002270 

0  vi-2280 

0.92  290 
0w23tr0 
05:2310 
U9232C 
052230 
0  i.  234 1) 
0»235& 
0-2^60 
.1)02370 

(j-2ii8i! 

oO2o90 
uu24l)G 
i)  u24i0 


i- 


I  ’ 


? , 
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EXP!  =  EXP  (~R*XL*CS) 

EXP2  =  EXP(  R*XL*CS> 

BETA  =  X  L‘*R<  SS 
GO  TO  CIO,  20),  L 
lu  THETA  =  Y*T.T3  -  PI/2.0 
SINN  =  3 IK (THETA  -  BETA) 

SINK  =  SIN (THET  A  4  BETA) 

COSH  =  COS ( THET  A  -  BET  A) 

COSO  s  COS (THET A  4  BeT  A) 

FUNK3  =  RQ,  /  (TERM2SQ  '  YSC)  * (EXPT* (Y*COSH  4  T£FH2*SINH) 
1  -  EXP2  *  (Y’COSP  4  TERM2*SINP)  ) 

RETURN 

Y*TT3  -  PI 


<Y  4  OM)**2)*(TERM2SCi  4  (Y  - 
:RH230)'*C0SH  4  2(j4(s*Y*AL*SINH> 


p 

SIHH  =  3 IN(THET A  -  BETA) 

K 

S1NP  =  SIN (THET A  4  BETA) 

i5 

COSM  =  COS (THET A  -  BETA) 

COSP  =  SOS (THET A  4  BETA) 

'H 

FUNX3  =  OfiYPQ.  /  ( (TERM2SS  4 

fi 

1  (EXPT*IRYSC  -  OH*  * 2  -  Ti 

4? 

2  EXS2H:(YSa  -  OH **  2  -  T 

RETURN 

END 


FUNP(PI) 

♦  H  /  3.0 


Eliu 

Q5*5***55*****:*.4.1 ,4**5  **Tfr**********ir***4*.5  *4  5*  *****.*.»*-*  4  *  ******* 

4  4  5**4  **♦*  *-* *****  ***55*5  ****  ******* 

FUNCTION  FUNR(T) 

COMMON  AL,  OM,  XL,  L,  PT,  ZO 
COHSON/TT  4/TT4 
TT4  =  T 
H  =  -PI  /  2.0 
N  -  1 
7HO  =  0i«6 
|<OUR  =  FUNP(H) 

ENDS  =  FUNP(O.O)  4 
SUHO  =  (ENOS  4  4. 3*F3UR) 

10  H  *  H  (  Z*  0 
N  =  2  *  N 
TWO  =  TWO  4  FOUR 
FOUR  =  0.0 

Y  =  H 
1=0 
1=14  i 

FOUR  =  FOUR  4  FUNP (Y> 

Y  =  Y  4  H  4  H 
IF  M  »L-T*»  M>  30  TO  20 
FUNR  =  (ENOS  4  2.0*TWO  4 
IF(A8S(SUHO  -  EUNR)  .LT.  1.0 E-5) 

SUMO  =  FUNR 
GO  TO  10 

END  ... 

**»  *****************  5*  5*5*54***55*4**4**4444* 

-JIM!  544  45-4  *  5  444*  *4  *  *4  4  54*  5  44 45  **4  55*  4*5***454  5 55**  44 44*45**5** 

FUNCTION  FUNP(Y) 

COMMON  At,  OH,  XL, 

COMMON/ TT4/TT4 
COHMON/D/O 
SINE  =  >IN(Y) 

COSINE  =  COS  (Y) 

R02  =  SORT  (  (1.0  4  0*C0SINE)  «- 
R03  =  SiRT(  (ZO  -  0*Al*SINE>**2  + 

R04  =  SIRTT  (20  4  0*AL*$INE)**2  4 


20 


4.0 ’FOUR)  *  H  /  3.0 
RFTURN 


L,  BI»  ZO 


002420 
00243Q 
002440 
002450 
002460 
00247-0 
0M2480 
012490 
Ou‘250  0 
UG2510  ; 
0u2520 
002530 
002540 
On 2550 
002560 
002570 
u02580 
0M>**2)  )  *01*2390 
002600 

)  002610 
002620 
002630 
*****44**002640 
4**44****002650 
002669 
00267-0 

O02680 

002690 
U027O0 
002710 
002720 
0U2730 
002740 
0  02750 
O02760 
Gu  2-7-7- P 
0  0278  0 
002790 
0u2  800 
QgESlO 
002820 
0>j2830 
00280a 
00285.0 
002860 
002870 
002880 
00283P 
0u29u0 
*5***5***0O23l0 
********  *1,02920 
002930 
002940 
002950 


<D*SIH£)^*2J 
A L* *2* (6.5  4 
A L* *24(0.5  4 


D4C03INE)* 
D*CQ5IN£>  * 


*2) 

*2) 
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Gl  =  COSINE 

G2  =  (l.U  4-  D*S1*  7  R02 
G3  =  AL  *  (u  .5  4  0*61)  /  R03 
G4  =  AL  *  (0.5  *  0*G1)  7  R04 
RC1  =  SaRTC  (1.0  +  Gl)  /  2.3) 

RC2  -  32RT  (  (1.0  +  G2)  7  2.0) 

RC3  =■  S3RTr(  tl.O  +  G3)  /  2.0) 

RC4  -  SQRT  (  (i.dO  4  G4>  7  2.0) 

RSI  =  SORT (  (1.0  -  Gl)  7  2.0) 

RS2  =  SORT (  (1.0  -  G2)  7  2.0) 

R33  =  SORT (  (l.U  -  G3)  /  2.0) 

RS3  =  —  RS3 

RS4  -  SORT (  (1.0  -  G4)  /  2.0) 

ACT  =  R  Ci  *  RC  2*>RC  3  *  RC  4  +  RCl»RC4»RS2*lRS3  -  R32*RC4*RSV*RS3  + 

1  R03*RC4*R5i*RS2  -  RSi» RS4*RC2*RC3  r  RS2*RS4»PC1^RC3  - 

2  R33*RS4»RC1*RC2  -  RSi*RS2*RS3“RS4 

AST  =  RS1*RC2*RC3*RC4  -  RS2*RCi*RC3*RC4  4-  RS3*RCl*RC2*R04  4- 

1  RS1*RS2*RS3*RC4  *■  RC1*RC2*RC3*RS‘4  f  RC1*RS2*RS3*RS4  - 

2  R02?RSi*RS3*RS4  *  RCS'RSl'RSa *RS4 

ARR  -  SORT (0  *  R03  *  R04  7  R02)  ♦  XL 
ARC  =  ARR  *  ACT 

ARS  =  ARR  *  AST 

FUNP  =  (-'l-i.p/PI)  *£XP(D* AL*TT4  * COSINE  -  ARC) *CpS <0*AL*TT4*SINE 

RETURN 

ENO 


003450 

003060 

0U3470 

003080 

003U30 

0U3100 

0u3110 

0U3120 

0U3130 

003140 

003150 

0U3160 

003170 

003180 

003200 

003210 

-ARS) 

003240 

003250 
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